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Abstract 

A theory is presented for modeling the evolution of rotor wakes as a function of axial distance 
in swirling mean flows. The theory, which extends an earlier work to include arbitrary radial dis- 
tributions of mean swirl, indicates that swirl can significantly alter the wake structure of the rotor 
especially at large downstream distances (i.e., for moderate to large rotor-stator spacings). Using 
measured wakes of a representative scale model fan stage to define the mean swirl and initial wake 
perturbations, the theory is used to predict the subsequent evolution of the wakes. The results in- 
dicate the sensitivitv of the wake evolution to the intial profile and the need to have complete and 
consistent initial definition of both velocity and pressure perturbations. 


Introduction 

The advent of high bypass ratio turbofans has 
made rotor-stator interaction noise one of the ma- 
jor sources of aircraft engine noise. This is true 
whether one is concerned with the tone compo- 
nent or the broadband component of noise spec- 
trum. Reduction of the rotor-stator interaction 
noise, therefore, has become an integral part of 
many engine noise reduction strategies. Naturally, 
the success of a reduction technique depends crit- 
ically on how well the underlying mechanism of 
the noise generation is understood. It is gener- 
ally agreed that the principal mechanism of rotor- 
stator interaction is the impingement of rotor flow 
unsteady perturbations on the downstream stator 
vanes. In the case of tone noise, the perturbations 
are the coherent (periodic) wakes of the rotor, and, 
in the case of broadband noise, they include flow 
turbulence ingested and/or produced by the fan. 

Over the past few decades a number of an- 
alytical and numerical methods have been devel- 
oped to predict the rotor-stator interaction noise 
by modeling the rotor wakes as specified upstream 
vortical perturbations and the stator as a 2D rec- 
tilinear, or a 3D annular, cascade of airfoils. For 
unloaded cascades, the common approach is to 
use linear unsteady aerodynamic theory wherein 


flow nonuniformities are treated as imposed up- 
stream disturbances purely convected by a uni- 
form mean flow but otherwise not affected by it 
(i.e., the frozen gust assumption). The noise field 
is then calculated as the farfield potential response 
of the cascade to the impinging vortical distur- 
bances. For loaded cascades, methods based on 
the rapid distortion theory are used to account for 
the effect of mean flow distortion of the upstream 
disturbances. Both the wake evolution and ra- 
diated sound are then determined as part of the 
solution [1]. However, to date, only 2D versions 
of the theory have been implemented in practical 
noise computation schemes. 

In both types of approaches the upstream 
and downstream mean flows are taken to be uni- 
form. In this case, the modal composition of the 
duct acoustics, which depends only on the mean 
flow and duct geometry, can be readily deter- 
mined. The magnitude of the acoustic modes, 
then, depends on the unsteady aerodynamic re- 
sponse of the cascade. For annular geometries, 
strip theory is often used to compute the cas- 
cade unsteady aerodynamics by treating the 3D 
response as an aggregate of 2D reponses at each 
radius strip (see, for example, [2]). 

However, despite the obvious simplicity of 
the mathematical development when the mean 
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flow is assumed uniform, there are definite short- 
comings in these types of approaches. This is 
mainly because the flow downstream of moder- 
ate to highly loaded fans show significant swirling 
motion. In fact, the swirl component of veloc- 
ity, which depends on the radial distribution of 
the blade loading, tip speed and rotor cascade 
stagger, may be comparable to the axial compo- 
nent of velocity in these flows. Wakes and sec- 
ondary flow nonuniformities, therefore, evolve in 
a swirling flowfield. The mean swirl produces cen- 
trifugal and Coriolis forces which affect the motion 
and structure of the perturbations present in the 
flow and couple the vortical, entropic and acoustic 
modes [3]. The magnitude and phase of the un- 
steady perturbations evolving in swirling flow are, 
therefore, constantly changing. 

Essentially, mean swirl has three principal 
effects on the rotor-stator interaction noise. First, 
it affects the evolution of the rotor wakes down- 
stream of the fan and, thus, changes the upwash 
velocity induced on the outlet guide vanes (OGV) 
[4]. Second, it modifies the aerodynamic and 
aeroacoustic response of the OGV [5]. Third, it 
changes the structure of the acoustic duct modes 
and their cut-off criterion [6, 7, 8], 

In trying to address the first issue, Golubev 
and Atassi [4] considered a mean swirl consisting 
of free- vortex and rigid body rotation components. 
Their results indicate that the influence of swirl on 
the wake evolution depends on the particular ra- 
dial distribution of swirl considered. For example, 
a free-vortex swirl distribution leads to algebraic 
growth of the wake perturbations. On the other 
hand, a general swirl distribution may lead to ax- 
ial modulation, exponential growth or decay of the 
wake perturbations. There is some experimental 
evidence for growth of wake perturbation under 
certain conditions as shown by the data of Laksh- 
minarayana [9] who measured an increase in the 
radial component of wake perturbation as a func- 
tion of downstream distance. It should be empha- 
sized that the growth or decay effect being dis- 
cussed here is separate from the viscous decay ef- 
fect which is always present in a real flow. The 
overall evolution of wake perturbations will de- 
pend on both these effects. However, in what fol- 
lows we shall ignore the viscous decay effect with 
the understanding that one can account for it, at 
least approximately, using wake decay correlations 


that abound in the technical literature. 

The objective of this paper is to examine 
the effects of mean flow swirl on the structure of 
the rotor wakes as they convect downstream. To 
this end, we first extend the analysis by Golubev 
and Atassi [4] to accommodate a general radial 
distribution of mean swirl present a model for the 
wake evolution. Second, flowfield measurements 
obtained at NASA Glenn Research Center for a 
typical fan stage scale model will be used to de- 
fine the mean swirl and initial wake perturbations 
for use in the theory. Third, data-theory com- 
parison of the wake evolution will be presented. 
Finally, the paper will conclude with a summary 
of important results. 

Theoretical Model for the Wake 
Evolution 

As stated earlier, the primary objective of the 
present paper is to account for the effect of mean 
flow swirl on the evolution of wake perturbations 
downstream of the fan. We are particularly inter- 
ested in describing the modification of the struc- 
ture of the perturbation flow as it evolves over a 
distance on the order of the mean duct radius r m . 
It is therefore plausible to assume that the fluid is 
inviscid and non-heat conducting. In addition, we 
suppose that the flow velocity and pressure can be 
written as 

V{x,r,9,t) = U(r) + u{x,r,6,t), (1) 

p{x,r,d,t) = po{r) + p'(s,r,0,t), (2) 

where x and r are the axial and radial coordi- 
nates, respectively, 6 is the circumferential an- 
gle, and t is the time. We further assume that 
\u(x,r,0,t)\ « \U(r)\ and |p'(x,r,M)| « |po(r)|, 
and that the mean flow velocity has a zero radial 
component, 

U(r) = U x {r)e x + U s (r)e$, (3) 

where (e x , eg) are the unit vectors in the axial and 
tangential directions, respectively. 

It is advantageous to represent the pertur- 
bation velocity as the sum of a vortical part, u^, 
and a potential part, V</> [4], 


2 


u = V0 + U R . (4) 

Substituting equations (1), (2) and (4) into the 
Euler equations we get 


where o = ax—u}t-\-m6. The parameter A is signif- 
icant, since it determines the flow stability. When 
A is real, the flow is stable and, hence, wake per- 
turbations remain bounded. However, when A is 
imaginary, the flow can become unstable resulting 
in growing wake perturbations. 

Next, equation (6) is discretized and solved 
numerically assuming u R = u R on the right hand 


D 0 u 

Dt 


-ft 


+ {u R - V)<7 = -(V xO)x V</>, (5) 


( Rl f Ly . (p 0 V)^ 4> = —V • ( p 0 u R ), (6) side. This provides the first iterative solution for 

\Dt \Cq Dt J p 0 ) Po S' 


where the operator and the pressure are given 

by 


Do 

Dt 


= 4 + ^v, 


dt 

P 1 = 


Dq4^ 

’ Dt' 


(7) 

(8) 


A particular solution for the vortical veloc- 
ity can then be calculated by solving, 


D °“” + <n R ■ VW7 = - (V x U) x V(f>. ( 12 ) 

Dt ' v ’ V ) 


Solution of equations (5) and (6) requires an ap- 
propriate set of initial and boundary conditions. 
The velocity components are specified at the in- 
let of the computational domain together with an 
impermeability condition on the walls of the duct. 
The pressure is specified at both the inlet and out- 
let of the domain. 

Following [4], we write 


The vortical velocity is then given by u R — u R + 

Equations (5) and (6) are solved iteratively 
until appropriate convergence criteria are met for 
both V0 and u R . 

Boundary and Initial Conditions 


roc 00 _ 

{u R ,<j)}{x,r,e,t) = / ]T {A m , ipm}{ x i r ) (9) 

00 

xe i(-u .t+mO+ax)^ 

where a is defined by 

— (— u)t + m9 + ax) = 0. (10) 

We begin by solving the homogeneous part 
of equation (5) which can be explicitly integrated 
using the initial data specified at x = 0. The so- 
lution, denoted as u R , is given by 


The boundary and initial conditions necessary for 
carrying out the solution procedure are as follows: 


At the inlet, x = 0, we specify 

1. u% = Uinitial, i-e., the initial data given 
from the measurment. 

9 dl/j — ~Pjnitial 
' &Z PoUx 

3. u R = -V0. 


At the exit x = Z, we take 
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• On the tip and hub radii, we take 

1. u?+^=0, 


where <p = tp{x,r)e ia . Often, the pressure is not 
measured at the inlet of the domain. As a re- 
sult, the initial conditions for the pressure must 
be imposed artificially by assuming either p = 0 

»'!=«■ 
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Experimental Flowfield Measure- 
ment 

The experimental wake data used in this paper 
were obtained at the NASA Glenn 9-foot x 15- 
foot Low Speed Wind Tunnel during testing of a 
high bypass ratio subsonic fan stage model which 
was designed and built by the Allison Engine Com- 
pany. Figure (1) shows a photograph of this model 
installed in the tunnel. The fan has a diameter 
of 22 inches with 18 blades, a hub-to-tip radius 
ratio of a = 0.5 and a design rotational speed of 
10,400 EPM. A Laser Doppler Velocimeter (LDV) 
was used to map the velocity field downstream of 
the fan at three axial positions (locations 1, 2 and 
3 shown in figure (2) corresponding to 1.2, 2.4, 
3.4 inches downstream of the fan trailing edge, 
each at three different corrected fan speeds corre- 
sponding to the approach (5,200 RPM), cutback 
(8,900 RPM) and sideline (10,400 RPM) condi- 
tions. The LDV data was acquired with the sta- 
tors mounted in the aft position as shown in fig- 
ure (2) and with a bellmouth inlet installed on the 
model (not shown). Only data acquired at the ap- 
proach condition are used in the analysis and the 
comparison is shown in this paper. 

A one-component LDV system was used 
during separate wind tunnel runs to measure the 
axial and tangential components of velocity down- 
stream of the fan. The rotor tangential position 
corresponding to each velocity measurement was 
determined using a once-per-rev signal fed into a 
shaft angle encoder. The encoder divided each ro- 
tor revolution into 900 angular bins. For the 18 
bladed fan this corresponds to 50 bins per blade 
passage. Each time a velocity measurement was 
made the rotor position was determined by sam- 
pling the output of the angle encoder. Typically 
some 42,000 velocity measurements would be ob- 
tained at each radial-axial location from which 
passage averaged (periodic) velocity profiles could 
be constructed. The reduction process from raw 
data to average profiles is illustrated in figure (3). 

Consider the raw data for axial component 
of the velocity at station 1 and at 35% span(radial) 
location shown in figure (3a). Note that, the ab- 
scissa in this figure, labeled the rotor position, can 
be equivalently treated as time or tangential loca- 
tion. The first step in the reduction process is to 


compute the once-per-rev average of the velocities 
occurring within each of the 900 bins. Figure (3b) 
shows the resulting once-per-rev average velocity 
distribution. The next step is to average the once- 
per-rev average velocity distributions that include 
all of the 18 blade passages into one representa- 
tive passage. This process smooths out blade-to- 
blade variations in the flowfield that might have 
occured due to manufacturing differences in the 
blade. The result, which may be called a passage- 
averaged representation, is shown in figure (3c). 
Note that, at this point both the incoherent part 
of the velocity field (i.e., turbulence) and any co- 
herent part that is locked to one rotor revolution 
have been taken out of the velocity data. 

Contour plots of the passage-averaged axial 
and tangential velocities are shown in figure (4). 
The view is from downstream looking upstream 
with the fan rotating clockwise. The passage con- 
tours have been repeated to illustrate the periodic 
transition from one passage to another. Low ve- 
locity areas appear as dark regions while high ve- 
locity areas appear as bright regions. These con- 
tours provide an overall view of the variations of 
the rotor flowfield in axial, radial and tangential 
directions. In particular, at a given axial station 
and a fixed radial position, the tangential vari- 
ations clearly show wakes of the fan blades. In 
order to use the theory developed in the previous 
section, the passage-averaged flowfield needs to be 
further separated into a mean component and a 
perturbation component. The mean part at each 
axial-radial position is obtained by a circumferen- 
tial average of velocity components. The pertur- 
bation part is then simply what remains once the 
mean is subtracted out. Finally, in accordance 
with the normal mode representation discussed 
in the previous section, the perturbation part is 
Fourier-decomposed into its constituent harmon- 
ics. 

Data-Theory Comparisons 

Before proceeding with the comparisons, it is 
worth analyzing the nature of the measured mean 
flow variations. The mean axial and swirl compo- 
nents are shown in figure (5a). As noted earlier, 
only results for the approach condition (i.e. 5,200 
RPM) will be presented in this paper. Both the 
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axial and swirl components show marked varia- 
tions with the radius as well as the axial locations 
downstream of the fan. This clearly demonstrates 
the inadequecies of the uniform mean flow assump- 
tion often used in modeling the wake evolution 
behind a fan. The observed variations are mainly 
due to flowpath contraction (see figure (2)). The 
viscous effects also contribute to the variations, 
but are likely to have only a small influence for 
the axial distances considered here. The distribu- 
tion of A is shown in figure (5b), which indicates 
that the wake is a stable bounded perturbation. 
Equation (11) shows that the wake defect magni- 
tude is modulated with an envelope of wavelength 
_A_ 

IV 

The first Fourier harmonic (i.e., m = 18) of 
the measured axial and tangential perturbations 
sure shown in figure (6). For these components, 
there is a significant drop in the amplitude of the 
harmonic between locations 1 and 2 everywhere 
along the radius. The corresponding change be- 
tween locations 2 and 3 is much smaller except 
near the tip. Viscous decay notwithstanding, the 
disproportionate changes between locations 1 and 
2 as compared with the changes between locations 
2 and 3, is likely caused by the proximity of loca- 
tion 1 to the rotor trailing edge which is only 30% 
of the chord upstream. The change in the struc- 
ture of the harmonic between locations 1 and 2 
indicates that there are significant flow variations 
in the tangential as well as the axial directions. A 
similar effect has been predicted in the theoretical 
analysis of reference [10] and the numerical predic- 
tions of reference [11]. It is important to note the 
radial redistribution of the tangential component 
as it evolves downstream (the harmonic appears 
to decay for r > r m and amplify for r < r m where 
Vffi = ( rhub + r tip)l 2). This effect will be shown to 
be caused by the mean swirl. 

Using the mean and perturbation quantities 
obtained from the measurements as initial data, 
the theory described in this paper is used to pre- 
dict the evolution of individual harmonic compo- 
nents of the perturbations as a function of the 
downstream distance. The results for the evolu- 
tion of the first harmonic (m = 18) are summa- 
rized in figure (7). Slices of these 3D plots for 
the axial, tangential and radial components are 
shown in figure (8). The first three downstream 
locations shown correspond to the measurement 


locations 1, 2 and 3, while the 4th location is fur- 
ther downstream at 7.4 inches downstream of the 
fan trailing edge. A comparison between the the- 
oretical evolution and the measured evolution for 
the first harmonic axial and tangential perturba- 
tion components is shown in figure (9). To obtain 
the complete flow field, the individual harmonics 
are reassembled as shown in figure (10). The com- 
plete flowfield may be used to compute the acous- 
tics [5].The descrepencies between the measured 
and predicted evolutions for the m = 18 harmonic 
component (see figure (9) may be explained as fol- 
lows. 

First, since only axial and tangential vel- 
coity components were measured, the initial distri- 
butions of the radial velocity and pressure pertur- 
bations must be artificially imposed on the goven- 
ing equations. Therefore, any errors in this speci- 
fication are likely to persists (or even grow) as the 
perturbations evolve dowsntream of the location 
1. To ensure a proper and consistent specification 
of initial data, therefore, it is important to have 
knowledge of all three velocity component pertur- 
bations as well as the pressure perturbation. More 
will be said about this later. 

Second, as was mentioned before, the influ- 
ence of the unsteady potential field of the rotor 
on the locations near the fan trailing edge is not 
included in the theory. Hence, any such influence 
is missing from the evolution equations. There 
are experimental evidence (see reference [9]) that 
a significant decay can occur in the amplitude of 
the wake for x/c < 0.6 due to the fan potential 
flow interference effects. 

Returning to the data-theory comparisons 
shown in figure (9), note that swirl mainly affects 
the disturbances near the tip section where the 
mean flow has large gradients. This can also be 
seen by observing the variation of A shown in fig- 
ure (5b). The predicted tangential perturbations 
are seen to have the same characteritics as those 
that were experimentally measured with a redis- 
tribution occuring around r ~ r m , as seen in the 
inset of figure (8b). 

The predicted radial component of pertur- 
bation is seen to amplify essentially monotonically 
as a function of downstream distance near the tip. 
Since, no measurements were taken for the radial 
component of the velocity, it is not possible to as- 
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sess this effect by comparison with data. There is, 
however, some evidence that the radial perturba- 
tion component is at least as large as those of the 
other two perturbation components lending some 
support to the predicted results [13, 14]. 

To demonstrate the importance of having a 
consistent set of initial data, calculations were car- 
ried out using an assumed initial profile of the ra- 
dial perturbation at x = 0 (i.e., locatin 1). The 
profile used in this calculation corresponds to the 
one obtained at location 4 in the previous calcu- 
lation. The resulting evolutions show significant 
difference compared with those for the zero initial 
radial perturbation case as shown in figure (11), 
although the effect of the swirl at the tip section 
is similar between the two cases. It is remark- 
able that the radial perturbations do not amplify 
anymore in contrast to the results shown in figure 

( 4 . 

We now examine the effect of modifying 
both the homogeneous pressure and radial velocity 
conditions at the inlet. The results show signifi- 
cant differences with the data and the other two 
cases. This only serves to emphasize the impor- 
tance of a consistent and compatible set of initial 
data containing all three velocity components and 
pressure distributions at x = 0. 

The specification of pressure perturbations 
at the inlet raises an important question regarding 
the boundary conditions imposed on <j) for the nu- 
merical solution. One might specify the pressure 
perturbations (as was done above), or one might 
choose to specify a derivative conditions such as 
= 0. The effect of the latter condition, which 
may be more realistic from a physical point of 
view, however, is rather small. This is expected 
because the vortical modes have a very small pres- 
sure content in them as pointed out in [4, 8] since 
these perturbations are nearly convected. 

At the exit of the computational domain, 
different numerical conditions were tested. 



The first condition is a common condition used in 
typical CFD calculations. The second condition is 


a weaker form of the first condition and is tested 
to check the importance of the term in the 
third condition. The third condition is based on 
the physical argument that the pressure pertur- 
bation is modulated by a wavelike structure with 
A- representing the associated wavenumber. The 
results of the analysis indicate that: 

1. The decay of the modes at the tip section 
was less for the first condition than the third. 

2. The choice of the second condition showed 
nearly the same results as those in the first 
with very little noticeable differences. 

3. The third condition yielded a smooth pres- 
sure profile at the exit. 

4. For all cases, the solutions showed similar 
basic trends with only small difference in de- 
tails. This lends support to the argument 
that the flow is a nearly convected flow as 
shown in [4, 8]. 

5. It may be concluded that the choice of the 
most physical condition is best for this type 
of evolution problem as was also suggested 
in [4], 

Conclusions 

The theory developed by Golubev and Atassi [4] 
was generalized to account for arbitrary radial dis- 
tributions of mean swirl in modeling the evolution 
of rotor wakes. Using experimentally meassured 
flow data as mean flow and initial condition in- 
puts, the theory was used to predict the evolution 
of the wake perturbations as a function of distance 
downstream of the fan. Comparisons with mea- 
sured wake evolutions on a harmonic basis indi- 
cate that swirl has a significant effect where large 
radial mean flow gradients occur. It appears that, 
for the most part, the radial and tangential pertur- 
bations exchange momentum as the perturbation 
field evolves downstream. 

The theory was unable to model the signifi- 
cant change between locations 1 and 2 as seen in 
the data. This may be attributed to the proximity 
of the measurement location 1 to the fan trailing 
edge where unsteady potential effects of the fan 
may be significant. 
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The effects of different inlet and outlet con- 
ditions on the evolution were also analyzed. It 
appears that the results are sensitive to the inlet 
conditions, but that outlet conditions play only a 
small role. This supports the normal mode analy- 
sis which indicate that the vortical perturbations 
have inly a small pressure content associated with 
them. 
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Figure 1: Allison Low Noise fan model in 9-foot x 15-foot NASA Glenn Low Speed Wind Tunnel 
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Figure 7: Numerical Simulation 
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Figure 8: 1 st Fourier mode evolution for the axial, swirl and radial velocities 



Figure 9: Comparison between the present theory and measured data 
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